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Abstract 

The constitutive behavior of aluminum alloy A356 in the as-cast condition 
has been characterized using compression tests performed over a wide range 
of deformation temperatures (30-500°C) and strain rates (~0.1-10 s" 1 ). This 
work is intended to support the development of process models for a wide 
range of conditions including those relevant to casting, forging and machin- 
ing. The flow stress behavior as a function of temperature and strain rate 
has been fit to a modified Johnson-Cook and extended Ludwik-Hollomon ex- 
pression. The data has also been assessed with both the strain-independent 
Kocks-Mecking and Zener-Hollomon frameworks. The predicted plastic flow 
stress for each expression are compared. The results indicate that the ex- 
tended Ludwik-Hollomon is best suited to describe small strain conditions 
(stage III hardening), while the Kocks-Mecking is best employed for large 
strain (stage IV). At elevated temperatures, it was found that the Zener- 
Hollomon model provides the best prediction of flow stress. 
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(3 strain rate function coefficient 

7 strain rate function coefficient 

e true strain 

e p total plastic strain 

e strain rate (s _1 ) 

Eq absolute strain rate (s _1 ) 

i a average strain rate (s^ 1 ) 

e* Johnson-Cook reference strain rate (s _1 ) 

T) Zener-Hollomon fit coefficient 

#i_2 Johnson-Cook strain rate adjustment coefficients 

strain hardening rate (MPa) 

Go initial strain hardening rate (MPa) 

Ki_2 Johnson-Cook thermal softening adjustment coefficients 

Ai_2 Johnson-Cook thermal softening adjustment coefficients 

A Johnson-Cook strain rate adjustment coefficient 

H temperature corrected shear modulus (MPa) 

/io absolute shear modulus (MPa) 

p Zener-Hollomon fit coefficient 

a recorded true stress (MPa) 

a predicted flow stress (MPa) 

a s saturation or Voce stress (MPa) 

ct s q absolute saturation stress (MPa) 

cr* Johnson-Cook reference stress (MPa) 
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<y y yield stress (MPa) 

$ Zener-Hollomon fit coefficient 

a,b,c extended Ludwik-Hollomon strain hardening function coefficients 

b burgers vector magnitude (m) 

C Johnson- Cook strain rate coefficient 

Cy,c y yield strength transition coefficients 

E temperature corrected modulus of elasticity (MPa) 

g normalized activation energy 

go initial normalized activation energy 

k Boltzmann constant (1.3806503 x 10~ 23 m 2 kg s" 2 K" 1 ) 

ki_ 3 extended Ludwik-Hollomon strength parameter function coefficients 

Kt Johnson-Cook thermal softening coefficient 

K(T) extended Ludwik-Hollomon strength parameter function 

M(T) extended Ludwik-Hollomon strain rate function 

N(T) extended Ludwik-Hollomon strain hardening function 

p dislocation interaction characteristic parameter 

Q Zener-Hollomon activation energy (J) 

q dislocation interaction characteristic parameter 

R gas constant (8.3144 J/mol- K) 

T instantaneous temperature (°C) 

T* Johnson-Cook reference temperature (°C) 

T a average temperature (°C) 

Th homologous temperature 
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T me \t melting temperature (°C) 

T t strain rate coefficient transition temperature (°C) 

Z Zener-Hollomon parameter 

1. Introduction 

Cast aluminum alloys continue to displace ferrous alloys in transport ap- 
plications where weight and part count reduction are important considera- 
tions. The manufacturing processes for commonly cast aluminum alloy A356 
(A17Si0.3Mg) components often includes machining operations prior to final 
heat treatment. Modelling and optimizing processes occurring before heat 
treatment requires a comprehensive characterization of the flow stress under a 
wide range of thermal and deformation conditions. Processing that occurs in 
this time frame affects the final component properties. This includes casting 
which involves low strains and strain rates at high temperatures, to forging 
which involves high strains and strain rates at intermediate temperatures. 

There have been many phenomenological and physically-based constitu- 
tive models that have been employed to predict plastic flow stress as a func- 
tion of material structure, temperature and strain rate. While a significant 
number of studies have been undertaken to study the deformation behav- 
ior of A356, all have focused on the material following various post-casting 
thermal treatments. The most common heat treatment for A356 is the T6 
heat treatment, which involves bringing the secondary phases into solution, 
quenching and then artificially ageing. As such, the deformation behaviour 
of A356 at room temperature has been well characterized in the T6 condi- 
tion [l|, 0, 0, 0] 5 an d at elevated temperatures in the solution-treated state 
The response of A356 during semi-solid forming [(J and rheocasting [?J 
have also been studied but the non-typical microstructure produced during 
these processes make this data of limited use for traditionally cast A356. 

In order to address the lack of deformation data for A356 at elevated 
temperatures and strain rates in the as-cast condition, a large number of 
isothermal compression tests were conducted in an effort to develop a compre- 
hensive constitutive equation. Prior to presenting the experimental results, a 
literature review of suitable constitutive expressions will be presented. Anal- 
ysis of the experimental data within each constitutive framework will follow 
with discussion of efficacy and performance. This will be accomplished by 
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comparing the experimental flow stresses to those predicted with the goal of 
identifying the best approach for a given thermomechanical state. 



2. Constitutive expression overview 



In general, the phenomenological deformation behavior of aluminum al- 
loys is dominated by strain hardening at low temperatures, which transitions 
to a time or strain-rate dependent (creep) response at elevated temperatures. 
Models of aluminum casting processes consider the material either from 
an elastic time-dependent creep standpoint, or as an elastic-plastic strain- 
dependent material. The former formulations considered the Sellars-Tagart 
j8[ or Garafalo relationships @, E3, Q, while the latter Q used 
extended Ludwik-Hollomon expressions. Constitutive behaviors of as-cast 
wrought aluminum AA1050, AA3104 and AA5182 have also been described 
with a simultaneous combination of time and strain dependence 15]. The 



Johnson-Cook [16( model, which is specifically structured to isolate strain- 
rate and temperature dependence, is another popular phenomenological ap- 
proach although it has not been employed for cast materials. Since there 
is often a co-dependance of temperature and strain-rate, many have modi- 
fied the original Johnson-Cook model to better describe specific materials at 



specific temperatures and strain-rate ranges [17|, |18|, |19| . While being quite 



effective in predicting flow stresses, the main drawback of using phenomeno- 
logical approaches is that they are only valid within the experimental strain 
range from which the constitutive expression has been derived. 

Physically-based models, such as the Zerilli- Armstrong [201 which con- 



siders grain size or the more general Kocks-Mecking |2l|, |22|, |23|, |24|, employ 
material parameters to predict flow stress. Grain size in cast Al-Si alloys has 
been shown to have little effect on str eng th compared to other microstruc- 



tural features such as dendrite spacing [25J, |26|, |27| . As a result, models which 



use grain size as their underlying basis are not well suited to material with 
dendritic microstructure. Based on its more fundamental derivation, the 
Kocks-Mecking approach has been selected as a candidate physically-based 
model for as-cast A356 in the current work. 



2.1. Kocks-Mecking 

The Kocks-Mecking hardening model is the most common physically- 
based approach to describe constitutive behavior. The underlying premise 
of the Kocks-Mecking hardening model is that the work hardening rate G = 
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da/ de approaches zero at some saturation stress, a s which is a function of 
strain rate and temperature. As a result, plotting 6//i versus a/a s creates 
a master curve that is accurate over a large range of hardening such that: 

I"© 

Typically, Eq. [T] is assumed to have the form of a linear Voce-type relation- 
ship: 

a 



e = e ( i - - ) (2) 



with an integrated form of: 



ex P — — ( 3 ) 



o y — a s 



The saturation stress, a s , is expressed as: 




(4) 



which is a function of a normalized activation energy term g is defined as 

kT e 

9 = ln — (5) 
fib 3 e 

where Eq is the minimum strain rate that converges all data to a single func- 
tion that relates cr s /fi to g. The physical basis of p and q (0 < p < 1 and 



1 < 9 < 2) in Eq. H] represent the shape of dislocation obstacle profiles 23 



which in turn decide values of <j s q and g^. These values are phenomenological 



in nature [24|, [28[, and have not been used to distinguish between discrete 
obstacle types such as precipitates or grain boundaries. The yield strength, 
cr y , can also be expressed as a function of g, which renders Eq. [3] a a-e 
relationship that is a function of strain rate and temperature. 

2.2. Johnson-Cook 

The Johnson-Cook model remains a widely adopted phenomenological 
model, potentially owing to simplicity and not requiring an extensive exper- 
imental regimen to develop. The heuristic basis of the Johnson-Cook model 



(3 



is such that the flow stress at elevated temperatures is determined by scaling 
the constitutive behavior at low temperatures and quasi-static strain rates. 
This approach consists of determining the product of three terms which is 
typically expressed as: 



a JC = a*(e)f(e)f(T) 




(6) 



The reference stress c*(e) in Eq. [6] is a relationship that describes the quasi- 
static flow stress at a low temperature as a function of strain, usually fol- 
lowing a standard Ludwik-Hollomon expression, cr*(e) = a* + K*e n * . The 
reference stress data is collected at reference strain rate e* (usually 1 s _1 ) at 
temperature T*. The second term describes the strain rate scaling based on 
a fitting coefficient C, and the third term is intended to capture the effects of 
temperature as a function of the homologous temperature, Th and a constant 



m (~ 1 for aluminum 16[) where: 



Th = t ( 7 ) 

J melt ~~ 1 

The main drawback of the standard Johnson-Cook model is that strain hard- 
ening is not explicitly affected by strain rate or temperature. The standard 
Johnson-Cook model has been modified by various functions of strain rate 



and temperature by other researchers to accommodate for this [29] . In this 
study, the Johnson-Cook framework has been utilized to describe the behav- 
ior of A356. 

2.3. Extended Ludwik-Hollomon 

Another phenomenological model frequently employed for cast aluminum 
alloys is the extended Ludwik-Hollomon expression. In order to capture 
the simultaneous evolution of both strain hardening and strain rate effects 



with temperature for aluminum, van Haaften |15| fit stress-strain data to an 



extended Ludwik-Hollomon expression of the form: 

,JV / .n .\M 



a L u = K{e + ey v {e° + e) M (8) 
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where K, N and M are constants corresponding to strength, strain hard- 
ening and strain-rate sensitivity, respectively. The constants e° and e° were 
reported by van Haaften to be numerically necessary in order to provide the 
correct yield strength at room temperature using the Simplex method for 



fitting. Expanding on the work of Phillion et al. |13l. Il4|. this expression has 



been modified such that K, M and TV are functions of temperature: 

M(T) 



<r LH = K(T)e 



N(T) 



( e V 1 ; 



where E\ is a normalization strain rate of 1 s -1 . The advantage of this 
phenomenological approach as opposed to the Johnson-Cook model is that 
N(T) and M(T) are both independent functions which when factored from 
the flow stress enable K(T) to be determined via regression. Compared to 
the single independent function that forms the basis of the Johnson Cook 
model, this is a more direct approach. 

2.4- Zener-Hollomon 

The Zener-Hollomon is an approach that has been widely used to char- 
acterize strain-independent, high-temperature flow stress. The flow stress is 
expressed as a function of the Zener-Hollomon parameter Z: 




where Z is expressed in terms of strain rate and temperature [30| : 



Z = f(e,T)=eexp[^-\ 



or stress by using the Sellars-Tagart |8[ Arrhenius-type equation: 

Z = $sinh(p(j) r ' (12) 

where Q is an activation energy, R is the gas constant, a is the strain- 
independent flow stress and fitting coefficients are p and rj. These latter 



three coefficients are occasionally reported to be solute-dependent [15] , how 



ever studies on steel have suggested that they are merely fitting parameters 



S 



and have no physical basis 3lj • The only parameter with a physical basis is 
the activation energy Q, which is an activation enthalpy. 

Owing to the wide range of constitutive expressions available, the preced- 
ing expressions will be applied to the experimental results described in the 
next section to determine the best assessment of the constitutive behavior of 
as-cast A356. 



3. Material and experimental methodology 

Al-Si-Mg alloys, such as A356, have an as-cast microstructure consist- 
ing of primary aluminum dendrites (cc-Al) which form during the initial 
solidification phase, surrounded by an Al-Si eutectic. As-cast microstruc- 
ture refinement, and a corresponding strength increase, is achieved primarily 
through decreasing solidification time or which decreases the secondary den- 
drite arm spacing (SDAS). Precipitation strengthening particles (Mg2Si) oc- 
cur throughout the microstructure, and other tertiary phases may be present 
due to impurities (Fe,C). The T6 heat treatment typically applied to this al- 
loy serves to modify the eutectic structure and refine/redistribute the Mg 2 Si 
particles. 

The material investigated in this study is a strontium-modified A356 cast 
from metal extracted from the melt supply of a North American aluminum 
alloy wheel manufacturer with the nominal chemical composition given in Ta- 



ble [H The melting temperature (T melt ) of this alloy is 612. 5°C [25]. Samples 
of this material were cast in a permanent steel mold with a water-cooled chill. 
Figure [1] shows an optical micrograph of the typical as-cast microstructure. 

Table 1: A356 composition in wt-% 



Element Si Mg Na Sr 

Range (wt-%) 6.5-7.5 0.25-0.4 -0.002 -0.005 



While the cooling rate was not measured during solidification, the SDAS 
and porosity were quantified with samples extracted at regular intervals from 
the casting. Quantification was conducted via optical microscopy and Clemex 
Vision PE software. Measurements were made at 12 locations representing 
a composite area greater than 85 mm 2 . The SDAS, defined as the average 
distance between the approximated centers of active secondary dendrite arms 
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Figure 1: As-cast A356 microstructure showing a-Al (bright), Al-Si eutectic 
(dark) phases. An example of a typical pore has been highlighted on the left 
and the high resolution of the microstructure on the right has been annotated 
to show the secondary spacing (SDAS). 



32] , was manually measured on both equiaxed and columnar dendrites. The 
average SDAS was found to be 39.9 \mi with a standard deviation of 7.0 
\xm. Both the average and standard deviation are based on 491 discrete 
arm measurements (Fig. [I]) distributed over the 12 locations; each loca- 
tion comprised at least 30 measurements. The measured gas-based porosity 
was widely dispersed with a mean percent area of 0.154%. This SDAS and 
porosity are representative of those observed in industrial castings from this 
manufacturer. 

Cylindrical compression test specimens, nominally measuring 10 mm in 
diameter by 15 mm in length, were extracted from the as-cast sample. De- 
formation tests were then conducted on a Gleebkfl 3500 thermomechanical 
simulator fitted with isothermal tungsten carbide anvils. A small amount of 
nickel-based lubricant was placed between the sample and the anvils. Dur- 
ing each test, the temperature was controlled with a type-K thermocouple 
mounted to the center of each specimen. Instantaneous diametral deforma- 
tion was measured with a linear voltage displacement transducer (LVDT) 



1 Gleeble is a trademark of Dynamic Systems, Inc., Poestenkill, NY. 
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actuated by quartz rods. This experimental arrangement is shown in Fig. |5J 
The procedure followed for each test was: 

1. The sample was first loaded between the platens with a pre-load of less 
than 0.5 kN. 

2. The hydraulic actuator was retracted by 2 mm to allow for thermal 
expansion. Pre-load on the sample was maintained via friction. 

3. The sample was Joule heated at a rate of 5°C s" 1 and then held at the 
target test temperature for 60 seconds. 

4. The sample was deformed at the prescribed deformation rate to target 
strains of 0.2, 0.4 and 0.55 based on the specimen's nominal length at 
ambient temperature. 

5. Based on the prescribed target strain rate, the data acquisition rate 
from the load cell, thermocouple and LVDT was greater than 200 sam- 
ples per unit strain during deformation. 




Figure 2: A cutaway depiction of the Gleeble experimental arrangement. The 
specimen with a thermocouple junction is positioned between the anvils and 
LVDT. 

At elevated strain rates, heat generated during the deformation was greater 
than could be removed by the anvils. This caused a slight increase (~ 10°C) 
in the average deformation temperature T a , and was more predominant at 
lower temperatures and higher strain rates. Although the pre-load at all tem- 
peratures was kept below 0.5 kN, creep occurred at elevated temperatures 



11 



during the hold stage prior to deformation. While the creep rate differs with 
temperature, the amount of creep was nearly identical from test to test as the 
pre-deformation heating cycle was the same. Differing amounts of thermal 
expansion caused by variations in specimen geometry render quantification 
of the amount of creep prior to deformation intractable. However, due to the 
low load and brief thermal history each sample experienced pre-deformation, 
this amount of creep is small and has been considered part of the thermal 
expansion included in the initial diameter, Dq. 

Uniform plastic deformation was assumed in all cases such that the instan- 
taneous true strain e was found from the instantaneous diameter D according 
to: 

, D 

e = -2 In — 

The instantaneous strain record was then differentiated to provide the in- 
stantaneous strain rate, e. The average strain rate, i a , calculated as the 
mean strain rate during deformation, was used to categorize each test. This 
is also true for the average temperature, T a , as well. All experimental data 
points of temperature, load, strain and strain rate have been incorporated in 
the following analysis. 



4. Experimental Results 

In total, 55 successful compression tests were carried out at a range of 
temperatures, strains and strain rates. Target strains and strain rates were 
calculated based on the nominal length of the specimen at ambient tempera- 
ture. Temperatures of 30°C, 100°C, and between 200 to 500°C in increments 
of 50°C and strain rates of 0.1, 1, 5 and 10 s _1 were targeted. Following the 
tests, the strain rates achieved were calculated based on the strain versus 
time information. The achieved strain rates based on each target strain rate 
were grouped in ranges of 0.06 to 0.11 s" 1 , 0.57 to 0.98 s -1 , 3.59 to 5.33 s" 1 
and 7.58 to 12.24 s" 1 , which are referred to further as R%, R2, R3 and R4, 
respectively. Flow curves are presented for each of these ranges of e a in Fig- 
ures [3J161 while all test conditions are summarized in Table |2j The variation 
in target strain rates and strains was attributed to thermal expansion and 
variation in specimen length. A number of tests at the highest target strain 
rates in each temperature range were not successful owing to non-uniform 
deformation and barreling. Tests where these conditions developed were de- 
fined as failed tests and their data is not shown nor used for subsequent 
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Figure 4: Characteristic measured stress-strain results for e £ i?2 
0.98 b" 1 . 
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Figure 5: Characteristic measured stress-strain results for e £ R 3 = 3.59 to 
5.33 s^ 1 . 
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Table 2: T a , e a , a y and e p attained for all tests 





W 


~ 30°C 






[b) ' 


- 100°C 




No. 


T 


4 


<7y 


£p 


No. 


T 






Ep 


1 


33.4 


0.06 


97 


0.13 


10 


103.7 


0.07 


95 


0.36 


2 


36.1 


0.07 


112 


0.37 


11 


102.4 


0.07 


104 


0.17 


3 


34.0 


0.07 


109 


0.18 


12 


104.3 


0.70 


120 


0.14 


4 


41.0 


0.57 


114 


0.42 


13 


112.4 


0.74 


102 


0.45 


5 


38.1 


0.71 


102 


0.19 


14 


106.9 


0.77 


106 


0.20 


6 


36.8 


0.71 


105 


0.19 


15 


104.7 


0.77 


117 


0.15 


7 


37.5 


0.72 


113 


0.19 


16 


106.6 


0.78 


111 


0.19 


8 


33.7 


0.74 


107 


0.14 


17 


121.5 


4.23 


121 


0.45 


9 


45.4 


3.95 


137 


0.35 


18 


126.9 


8.27 


117 


0.39 




(c) 


~ 200°C 






(d) 


- 250°C 




No. 


T a 




Oy 


£ P 


No. 


T a 


£o 


Oy 


Ep 


19 


201.3 


0.08 


85 


0.19 


25 


249.8 


0.08 


95 


0.42 


20 


199.5 


0.08 


84 


0.32 


26 


251.6 


0.08 


91 


0.20 


21 


204.9 


0.72 


84 


0.31 


27 


247.5 


0.08 


89 


0.33 


22 


203.7 


1.10 


99 


0.29 


28 


257.9 


0.70 


101 


0.44 


23 


188.7 


3.86 


105 


0.27 


29 


259.9 


4.39 


106 


0.40 


24 


201.6 


7.58 


111 


0.27 


30 


258.5 


8.24 


111 


0.30 










(e)~ 


300°C 














No. 






a y e p 









31 297.8 0.08 90 0.33 

32 300.5 0.09 85 0.21 

33 307.6 0.80 101 0.47 

34 310.2 3.59 108 0.34 
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Table 2: T a , e a , a y and e p attained for all tests (cont'd) 





(0 


- 350°C 






(g)' 


- 400°C 






No. 


T a 




Oy 


Sp 


No. 


T a 


£a 


<7y 




35 


350.1 


0.08 


60 


0.20 


40 


400.7 


0.09 


43 


0.22 


36 


345.3 


0.09 


65 


0.34 


41 


396.8 


0.09 


45 


0.38 


37 


355.0 


0.81 


71 


0.47 


42 


404.3 


0.83 


45 


0.50 


38 


360.6 


4.92 


84 


0.48 


43 


407.2 


4.28 


56 


0.48 


39 


351.9 


8.82 


81 


0.35 


44 


407.2 


8.59 


71 


0.35 












45 


408.7 


9.17 


59 


0.37 












46 


404.4 


9.68 


61 


0.15 




(h). 


- 450°C 








(i)' 


- 500°C 






No. 


T 




er. y 


£ P 


No. 




£a 


(Jy 


£ P 


47 


AAA. A 


0.10 


27 


0.41 


51 


494.9 


0.10 


17 


0.40 


48 


452.5 


0.88 


33 


0.52 


52 


502.3 


0.98 


24 


0.54 


49 


458.6 


5.33 


42 


0.53 


53 


502.6 


4.97 


32 


0.55 


50 


445.5 


10.35 


50 


0.40 


54 


498.8 


11.21 


39 


0.46 



55 502.8 12.24 38 0.19 
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The yield stress for each test (shown in Fig. 3-6 as circle symbols on each 
curve) was found with a 0.2% offset method employing temperature-corrected 
shear and elasticity modulii 33[, /i and E, according to: 



Li = 2.54 x 10 4 (l + ^—^-) (13) 

V 2-Tmelt / 

E = 2li(1 + v) (14) 

based on pure aluminum, with T in Kelvin and v = 0.33. A relationship 
based on pure aluminum is appropriate as the temperature dependence of 



E for aluminum alloys is relatively insensitive to solute content |34|. The 
sparse nature of the experimental data in the elastic region of the flow stress 
coupled with the diametral strain measurement precludes using experimental 
data directly for temperature corrected modulii. 

The mechanical response of this as-cast material exhibits three temperature- 
dependent regimes of behaviour. Below 300°C, the flow stress of the mate- 
rial shows appreciable strain hardening at all strain rates and little strain 
rate sensitivity. Between 300-350°C, the strain hardening diminshes and the 
strain rate sensitivity increases owing to dynamic recovery. Above 350°C, 
there is little to no strain hardening and continued strain rate sensitivity. 
The strain rate dependence is also demonstrated by the trend in a y (Table 
12]). For the low temperature tests, a y remains relatively insensitive to strain 
rate as compared to tests at elevated temperatures. For the tests conducted 
at ~ 30 and ~ 100°C, a y for e a £ R\,R2 was found to be approximately 
108 MPa varying by ±10%. The estimated experimental error in the load 
transducer (±0.01kN) used in this study is two orders of magnitude lower 
than the variance observed in the measured yield stress, thus the majority of 
the variability in the measured yield stress is attributed to small variations in 
the SDAS and the coarse Mg 2 Si precipitates in the as-cast structure. Owing 
to the low area percent and the fact that the applied load will tend to close 
pores, porosity is not considered to influence the compressive strength of the 
material [Hf. 

For comparison, the low temperature a y observed in this st udy of as-cast 



material is half that of similar material in the T6 condition [35[ and twice 
that of the solutionized condition jEJ. The main cause of this difference is 
surmised to be the size and distribution of the Mg 2 Si precipitates. In the 
solutionized condition, the precipitates have dissolved forming a supersat- 
urated solid solution. This results in an ineffective medium for dislocation 
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storage and a substantively smaller a y compared to the as-cast condition. 
The artificial ageing step in the T6 process serves to form smaller, more 
evenly distributed Mg2Si precipitates as compared to the coarse variety in 
the as-cast condition. This in turn serves as a more effective dislocation 
storage mechanism, resulting in a higher a y . 

For tests conducted at the lowest strain rate, R±, there is evidence of slight 
strain-softening at temperatures near 300°C (Fig. which diminishes at 
higher temperatures for tests in the same strain rate range. This behavior was 
also observed in the 300°C test for R 2 (Fig. [5]). These results are evidence of 
the dynamic phenomena known to occur for this alloy system at these temper- 



atures [38[ such as recovery and strain ageing. As 300°C has been identified 
as the start temperature for the transition from strain-dependence, the mate- 
rial behavior may be conservatively characterized such that below 350°C, the 
material is dislocation interaction dominated. Above this temperature, it is 
strain rate/diffusion dominated. As there is a temperature/strain-rate range 
where both dislocation interaction and diffusion processes are active, a con- 
stitutive relationship that traverses this regime will inherently have difficulty 



accurately describing the flow stress |29l . 138 



5. Constitutive equation development 

While the experimental results may be divided into two discernable regimes, 
from an analysis standpoint it is necessary to construct an equation or series 
of equations that successfully predicts the flow stress across both regimes. As 
outlined in Section [2j physically-based and phenomenological expressions will 
be applied to the experimental data to identify the most accurate. The ac- 
curacy of applying a physically-based expression to only the diffusion regime 
will also be presented. The analysis in the proceeding was accomplished 
through linear least squares fitting where possible and a Nelder-Mead tech- 
nique when non-linear fitting was necessary. The fitting procedure for each 
approach will be presented initially, after which the results of each model will 
be discussed. 

5.1. Kocks-Mecking 

As outlined in Section 12.14 the Kocks-Mecking description of constitutive 
behavior consists of relating a saturation stress a s to a normalized energy 
term, g. Values of a s were first identified from linear intercepts of — a data 
between O = and O = yu/20, which is defined as the transition from stage III 
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to stage IV |36|,|24|. The result of this process is demonstrated in Fig. [7] which 
shows Q/fi versus a/a s for strain rate ranges Ri and i? 4 . The 0//i versus 
a/a s demonstrates overall a continuous function satisfying Eq. [TJ While the 
large strain, high temperature data below O = yu/20 may be approximated 
by Eq. [21 the overall relationship is non-linear. Above this point, the data 
corresponding to small strains and low temperature represents a large initial 
work hardening rate. The initial work hardening rate 0o across all tests has 
been identified as 2875 MPa, as compared to annealed pure aluminum having 
6 = 1120-1720 MPa El ■ 
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Figure 7: Q/fi versus a/a s for e a G R\ and i? 4 at all temperatures. The 
horizontal line indicates = fi/ 20. 



Fig. [8] shows a plot of versus g(e a ,T a ) with p = 1/2, q = 2 and 
Eq = 10 7 s _1 . These constants are baseline values identified by Kocks and 
Mecking 24j. In Fig [71 the constants u s o/ and go may be determined from 
the y and x intercepts, respectively. A sensitivity analysis conducted on the 
values of p, q and Eq did not show appreciable improvement of the fit within 
the range defined by their physical basis. 

Fig. [8] also shows the relationship between a y and g. Coinciding with 
the two regimes identified in the previous section, there is an identifiable 
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Figure 8: Saturation and yield stresses versus normalized activation energy 
g, corrected for obstacle profiles by p = 1/2 and q = 2. 



threshold value of g where a y decreases; cr y /fi remains constant until this 
threshold is reached. For larger values of g beyond this transition, a y is 
taken to approach g . In order to characterize this transition, a bi-linear 
function was fit corresponding to data divided at g 1 ^ = 0.6: 



i/p 



g < (a + gl /q ] 



1/9 

^0 



where C y and c y are constants. The values of all constants are given in Table 

El 

As both a s and a y are able to be scaled as a function of temperature and 
strain rate, an expression for flow stress is possible through Eq. |3J However, 
this relationship is incapable of describing a large portion of the flow stress 
accurately, owing to the particularly high initial work hardening rate, and 
the non-linear strain-hardening rate at low temperatures and strains (Fig. 



[7]). Using a hyperbolic equation to describe the — a behavior 22] has been 



found to provide a better phenomenological representation of a versus e post 



21 



Table 3: Values of fitted coefficients in each material model 



Model 



Eq. 



Coefficients 
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O"s0 


//x = 0.0297 


#0 


= 1.1798 
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= 0.0656 
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= -0.1361 
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= 2850 MPa 
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El 
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EH 
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= 265.7 MPa 
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= 4024 MPa 






da 
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= 1.4125 
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= 1.4584 
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= 0.4143 


A 2 


= 4.5004 




loni 




= 0.312 


e 2 
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= 4.0072 
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= 346.3°C 
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= -5.8165 x 10" 4 


Extended 


EH 


6 = 


= 0.1962 


c = 


= -5.2734 x 10~ 3 


Ludwik- 


E2 


a = 


= 4.6489 x lO" 11 




= 3.5237 


Hollomon 


7 = 


= 0.2371 












= 8.1200 x 10" 4 


h 


= -1.1570 




E3 


^3 


= 407.7 






Zener- 
Hollomon 


El 


$ : 

77 = 


= 2.7209 x 10 19 
= 1.0265 


p = 


= 0.3700 


EH 


<? = 


= 2.3690 x 10 5 J 







yield for all conditions such that: 

(\ !/ 2 

Both in this expression and in Eq. [31 G>o is taken to be a static value, where a s 
and <j y are based on preselected dislocation parameters (p, q) for an absolute 
strain rate of Eq. This differs from other approaches taken to assessing these 
parameters for Al-Mg alloys, which involve regressively solving an objective 
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function [38] , which can potentially predict negative flow stresses under con- 
ditions resulting in low flow stresses. The other caveat of using a static value 
of Go is that the slight strain hardening followed by softening discernable 
at 300°C is not captured. However, this phenomena is only pronounced for 
a small range of thermomechanical states. Overall, this approach provides 
a reasonable description of the entire range of behavior observed without 
having to scale 0o- 

5.2. Johnson-Cook 



As described in Section 12. 2\ the Johnson-Cook constitutive behavior is 
based on scaling a reference flow stress expression to account for changes in 
strain rates and temperatures. In order to identify a reference flow stress 
expression, the test data was aggregated for samples tested at ~ 30° C (T* = 
34.3°C in Eq. [7]) and e 6 R\. Owing to the significant initial work harden- 
ing rate as described in the previous section, a single value of n* does not 
adequately describe the work hardening rate at all strains. Therefore, a hy- 
perbolic Voce-type saturation stress model of the same form as Eq. ITSl was 
substituted for a* (Eq. [6]): 

(\ 1/2 
{<? + (K) 2 - K) 2 ) exp (^J J (17) 

where a* and a* coincide with the average values for e a 6 R\ at ~ 30°C 
This approach expands on the linear Voce-type characterization that was 



previously applied by Alkarca et al. [39] in an amended Johnson- Cook type 



model for A356 in the T6 condition. The reference strain rate was taken to 
be the overall average strain rate of R\. 

The thermal softening function f(T) was derived by rearranging Eq. [6] 
to give: 

f(T) + f \T)C In (L) (18) 



Values for f(T) may be determined as the linear intercepts of oja* versus 
ln(e/e*). The values of f(T) were calculated from the reference tests for 
strains between 0.05 to 0.12 at increments of 0.01. This strain range was 
selected because it represents the post-yield response and terminates prior 
to the peak strain e v attained for all tests. 
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Comparing the standard form of f(T) = 1 — T^ 1 (Eq. E]) with the calcu- 
lated values of f(T) versus the homologous temperature T H indicates that 
m is greater than 1 at low temperatures where strain hardening is domi- 
nant, but is less than 1 at elevated temperatures where strain rate effects are 
dominant. As a result, a single coefficient of m does not adequately capture 
the thermal softening term. Therefore, the standard expression for f(T) has 
been expanded such that there are two terms to describe each of the ranges 
of behavior: 

f(T) = 1 - KiT* Al + k 2 T* x ' 2 (19) 

where K\ and Ai apply for low temperatures while k 2 and A2 describe elevated 
temperatures. The reference temperature T* was taken to be the minimum 
temperature found for all strains. A comparison between Eq. [TJJ]and f(T) = 
1 — is shown in Fig. IHH 

The strain rate coefficient C was then determined by dividing the slope 
term in Eq. ITSlby Eq. [T9j Plotting C versus Th (Fig. l9b|) shows that this 
term increases exponentially with temperature past a critical homologous 
temperature such that: 

C = 9 1 T£ + 9 2 (20) 

5.3. Extended Ludwik-Hollomon 

Outlined in Section [2731 and Eq. |9l the extended Ludwik Hollomon expres- 
sion for flow stress is coupled with temperature via strength, work hardening 
and strain rate parameters (K,N,M). Values of the work hardening param- 
eter, N, were found first via the slope of In a- In e curves. The strain-rate 
term, M, was taken to be the slope of In <r-ln e curves built from data at the 
same strains used in the Johnson-Cook analysis. The resulting values of M 
and N versus T a are shown in Fig. HOal These two data sets were then fit 
with Eq. [2T] and [22] and were used to extract values for K (Fig. UObl) . 

The N values plotted versus average temperature demonstrate that there 
is a roughly linear trend for all strain rates until approximately 350°C. Here, 
the diffusion effects identified in Section H] become significant. The values 
for e a G i?i are distinctly negative, while the rest of the strain rates result 
in iV values close to zero. Corresponding to regimes exhibiting strain-rate 
dependent phenomena, N(T) was characterized such that: 



JV(T)= i^i' (21) 
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_-/(T) = l-T™ 

f{T) =1-kiT£ +K2T£ 
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(a) Thermal softening behavior 
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lure 9: Johnson-Cook parameters plotted versus homologous temperature. 
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Figure 10: Extended Hollomon coefficients plotted versus temperature. 
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where T t is the transition temperature (~ 350°C) between the two regimes, 
with a, b and c constants fit using values of T a > 300°C and T a < 300°C 
respectively. 

As the physical basis for the trend in the M — T a data (Fig. I10a[) is an 
observed increase in strain rate sensitivity with temperature, a continuous 
function similar to that used for C (Fig. l9b|) has been employed for M(T): 

M(T) = aT? + 7 (22) 

The correlation of the strength coefficient with temperature, K(T), is con- 
tinuous and was found to be best described by: 

K{T) = k{T 2 + k 2 T + k 3 (23) 

5.4- Zener-Hollomon 

For the Zener-Hollomon analysis outlined in Section 12. 4[ it was assumed 
that the diffusion behavior at low strain rates is negligible. The Zener- 
Hollomon parameters were determined with T a > 350°C. The constants in 
Eq. [12] and [11] were fitted by starting with an initial value of p, which was 
then used to find 77, followed by Q and <3>. The value of p was then iterated 
upon to converge on a single value of Q at each strain between 0.1 and 0.15 
in increments of 0.01. Values of 77 were found as the reciprocal of the slope 
of ln(sinh(pcx)) versus ln(e) 0. By rearranging Eq. [TT1 and [T2l 

ln( S1 nh(p ff ))-^ = 4--!^ (24) 

77 7]R1 77 

Iterating over a range of p provides values of Z as a function of strain rate, 
temperature and stress (Fig. 11), converging at Q — 237 kJ/mol (Fig. 
[TT]) . In order to assess the sensitivity of the technique, the data from high- 
temperature tests (>~ 400°C) for select strain rates (e G -R2-4) was fit in 
the same manner. This data set was selected because it does not exhibit 
softening effects. The activation enthalpy, Q, for this data set was found 
to be 157 kJ/mol, which is considerably closer to the activation energy for 
self-diffusion of aluminum and diffusion of magnesium in aluminum, 131 and 
137 kJ/mol respectively. 
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Figure 11: Z = f(e,T) compared to Z = f(a). 



6. Comparison of constitutive expressions and discussion 

Fig. [12] shows examples of the measured flow stress compared to the 
predicted flow stress calculated by the different constitutive expressions. In 
these calculations, experimental data points for strain, strain rate and tem- 
perature are directly substituted into each respective constitutive expression 
for flow stress. All constitutive relationships predict the measured stresses 
reasonably well considering material strength variations (~ 10%) Figures fl2al 
and Il2bl demonstrates the characterization of the Kocks-Mecking, Johnson- 
Cook and Ludwik-Hollomon expressions at temperatures that exhibits strain 
hardening (100°C and 200°C), while Figures Il2cl and I12dl shows these ex- 
pressions as well as the Zener-Hollomon characterization for the strain rate 
sensitive regime (450°C and 500°C). Owing to the negligable strain rate sen- 
sitivity at 100°C, a single strain rate is presented corresponding to R\ in Fig. 
12 112al while the latter plots in Fig. 12 are representative of the lowest and 
highest experimental strain rate ranges, Ri and R4. 

At lower temperatures, all applicable relationships describe the experi- 
mental flow stress well. The Ludwik-Hollomon fits the data best for these 
conditions, fitting the data for both strain rates particularly well at 200° C 



28 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

e e 



(c) ~ 450°C (d) ~ 500°C 

Figure 12: Comparison of experimental versus predicted flow stress for ap- 
plicable models. 
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(Fig. I12bp . The Kocks-Mecking underestimates the low strain rate flow 
stress at both 100°C and 200°C, but is acceptable for the elevated strain rate 
at 200°C. The Johnson-Cook expression overestimates the flow stress in all 
cases, a trend that continues at elevated temperatures. 

At elevated temperatures, the Ludwik-Hollomon, Kocks-Mecking and 
Johnson-Cook expressions are fairly accurate at low strains, but overesti- 
mate the flow stress at elevated strain rates. The Zener-Hollomon expression 
captures the experimental data at all strain rates compared to the other 
three. Of the three expressions capable of capturing strain hardening, the 
Johnson-Cook is clearly incorrect even from a phenomenological standpoint: 
it demonstrates significant hardening even at 500°C (Fig. I12dp . 

Calculating the sum of the mean difference as well as the mean square er- 
ror (MSE) between the predicted flow stress using experimental data points 
and the experimental flow stress measured directly for all tests conducted 
provides a measure of each models performance. For the constitutive expres- 
sions applicable across both regimes, the Kocks-Mecking slightly underesti- 
mates the flow stresses by less than 1% with 9% MSE, the Johnson-Cook and 
Ludwik-Hollomon overestimate the flow stress by 8% and 2%, respectively. 
The overall MSE of the Johnson-Cook expression is 10%, however the MSE of 
the Ludwik-Hollomon expression is less than the Kocks-Mecking expression 
at 6%. The Zener-Hollomon model overestimates the applicable flow stress 
by 5% and has a MSE of 5% for all experimental data where it has been 
fit. In terms of specific strain rate and temperature ranges best described 
by each expression, the general trend for the strain-hardening conditions is 
that higher error is observed with increasing temperature, with peak error 
occurring in the transition regime (300-350°C). This is especially true for the 
Johnson-Cook expression, where the error reaches 25% for the 300-350°C 
range, while the others were approximately equal to their MSE values. 

As this material exhibits a range of temperature and strain rate dependent 
deformation response, selecting a single constitutive expression to capture 
the full range of behavior observed requires compromises in accuracy. Even 
though the Johnson-Cook model is widely used to generally characterize con- 
stitutive behavior, applied to as-cast A356, it is not an effective phenomeno- 
logical approach. This is despite the heavy modification applied in this study 
to improve its performance (Section 15. 2p . The extended Ludwik-Hollomon 
model provides additional flexibility in a phenomenological approach, how- 
ever requires a large number of fitting coefficients to approximate the flow 
stress effectively. The Zener-Hollomon technique describes the flow stress ef- 



30 



fectively, but is only applicable in the rate-dependent regime. The physically- 
based Kocks-Mecking is more efficient in terms of fitting coefficients, however 
it requires well characterized yield and saturation stresses. Furthermore, the 
accuracy of the Kocks-Mecking expression at low strains is dependent on the 
integrated form of Eq. [TJ in order to revert to consider strain dependence. 
The Kocks-Mecking expression in the integrated form may only describe the 
strain-softening behavior of the material by scaling the initial strain hard- 
ening rate according to temperature and strain rate. This requires a large 
number of tests to accurately assess, which may also be a limit when consid- 
ering inherent material variations. As an alternative to this physically-based 
model, the Ludwik-Hollomon approach offers a phenomenological description 
which has a greater degree of freedom to capture specific material behavior. 

7. Summary and Conclusions 

The constitutive behavior of as-cast A356 has been experimentally char- 
acterized through an extensive set of compression tests. The data was used 
to fit four phenomenological and physically-based constitutive expressions. 
The material in the as-cast form shows a diverse range of thermomechan- 
ical behavior characteristic of Al-Si-Mg alloys, transitioning from strain- 
hardening to being completely strain rate dependent. Below 350°C, the 
material was observed to strain harden, whereas above this temperature it 
becomes more strain-rate dependent with increased temperature as outlined 
in Section 3. This is also demonstrated in Section 5, most prominently by the 
work-hardening coefficient, N, plotted versus temperature (Fig. 10a) in the 
extended Ludwik model. Here, it is effectively at 300°C Furthermore, the 
inflection point of f(T) in the Johnson-Cook model (Fig. 9a) corresponds 
to 300-350°C Finally, g x l q = 0.06 in Fig. 8 also corresponds to this range of 
300-350°C, where a y approaches cr s . 

This behavior poses a challenge for developing a single constitutive ex- 
pression that accurately describes all experimental results across the tem- 
perature and strain rate ranges encompassed experimentally. As a result, 
each constitutive expression has differing degrees of accuracy and efficiency 
in predicting flow stress for particular ranges of temperature, strain rate and 
strain. Specifically, over the experimental data tested: 

• The extended Ludwik-Hollomon expression is the most flexible and 
therefore provides the best prediction across all temperatures and strain- 
rates. This expression overestimated the flow stress by an average 2% 



31 



with a MSE of 6%. This phenomenological approach necessitated a 
large number of fitting coefficients. 

• The Kocks-Mecking relationship on average slightly underestimated the 
flow stress by 1%, but with a larger MSE. Flow stresses for elevated 
strain and strain rate conditions showed better agreement with the 
relationship. This model is physically-based, and does not have as 
many fitted coefficients as the Ludwik-Hollomon. 

• The Johnson-Cook expression, did not describe the experimental flow 
stress effectively, even with heavy modification. Not only did the rela- 
tionship overpredict the flow stress, it exhibited strain-hardening at all 
temperatures. 

• For temperatures greater than 400°C where the material is completely 
strain-rate dependent, the Zener-Hollomon relationship shows a much 
better fit to experimental data at all strain rates as compared to the 
other models. However, the Zener-Hollomon expression is not valid at 
lower temperatures where strain hardening occurs. 

Based on this analysis, it has been found that the constitutive expression that 
most accurately describes as-cast A356 across all temperatures and strain 
rates is the extended Ludwik-Hollomon expression. 
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